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Abstract 

We study theoretically the temporal evolution and the spatial structure of the interface 
between two polymer melts involving three different species (A, A* and B). The first melt 
is composed of two different polymer species A and A* which are fairly indifferent to one 
another (Flory parameter xaa* — 0). The second melt is made of a pure polymer B which 
is strongly attracted to species A (xab < 0) but strongly repelled by species A* (xa*b > 0). 
We then show that, due to these contradictory tendencies, interesting properties arise during 
the evolution of the interface after the melts are put into contact: as diffusion proceeds, the 
interface structures into several adjacent "compartments" , or layers, of differing chemical 
compositions, and in addition, the central mixing layer grows in a very asymmetric fashion. 
Such unusual behaviour might lead to interesting mechanical properties, and demonstrates 
on a specific case the potential richness of multi-component polymer interfaces (as compared 
to conventional two-component interfaces) for various applications. 

1 Introduction and motivation 

The phenomena taking place at the interfaces and surfaces of polymeric systems are obviously 
of high importance to many practical or industrial situations. As a result, they have been the 
focus of extensive study both experimentally and theoretically in the recent decades. Much 
progress has been made in the understanding of these phenomena as well as in the development 
of techniques to characterize them (see e.g. the general references 1-4). 
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1.1 Interfaces between polymer melts 

One situation which has received much attention is that of two pieces of molten polymer put 
into (good) contact at a certain initial time: how does then the interface between these melts 
evolve, and how does the mixing between the melts occur, if any? or, more precisely, what is 
the final, equilibrium state reached by the system (i.e., given enough time, does it mix fully or 
only in a restricted region?), and with what kind of concentration profile? What is the dynamics 
leading to that final state, over which typical timescales? etc. It has been found, 1 ' 2 due to the 
specific physics of polymer macromolecules, that their mixing is generally quite different from 
what is known in more conventional diffusive systems, like molecular gases. 

Once the two polymer pieces are put into contact, the subsequent evolution will obviously 
generally depend on the various physico-chemical properties of the materials facing each other 
(nature of microscopic interactions between monomers, chain length, topological structure, . . . ). 
It has been recognized however that, at least from a conceptual standpoint, the two most crucial 
parameters determining the fate of the system are the sign and magnitude of the product xiV, 
where x is the Flory parameter, 5 related to the microscopic interactions between components, 
and N is the chain length. We will remind of the different possibilities for %iV in section 2. 

1.2 Multi-component interfaces 

In order to facilitate the identification and modelling of the fundamental processes at work 
within polymer/polymer interfaces, past studies have for the most part, if not exclusively, focused 
on interfaces involving two species (for instance, one on each side of the interface) and sometimes 
only one species (interfaces between two identical melts). Certainly, the span of real situations 
is much wider, with a vast range of multi- component interfaces, and one is thus entitled to ask 
what is the structure and dynamics of an interface between two polymer mixtures involving 
more than two components. 

As exemplified on the case studied in this article, such multi-component interfaces may have 
an interest of their own, as a means (for example) of obtaining structured interfaces at a sub- 
micron scale: when two polymer mixtures are put into contact, and provided the contrast in 
relative miscibilities of the different components is adequate, a "multiple" interface may form; 
that is to say the interface will partition into several adjacent layers, each of different chemical 
composition. A possible application of this could be the formation of multilayer polymer films. 
The mechanical properties presented by the resulting assembly might also prove interesting. 

Multi-species interfaces can also appear in the course of an interfacial reaction between two 
polymers: when the polymer pieces are put into contact, a reaction may start at the interface that 
delivers a (polymer) product whose miscibilities with each of the initial reactants are different; 
this three-species system is susceptible to lead to a multiple interface, whose dynamics is then 
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coupled to that of the chemical reaction. 

However, on purely combinatorial grounds, the diversity of all conceivable multi-component 
interfaces makes any general approach likely to be hopeless; thus delineating and focusing on a 
certain number of limiting situations may be useful. This is indeed the purpose of the present 
article: we will be considering one such "selected case" for a three-component interface, and will 
show that it exhibits a peculiar spatial structure and evolution through time, which we will fully 
analyze. 

The article is organized as follows: in section 2, we start with a reminder on the standard 
theoretical results concerning two-component interfaces, as these will be used as basic tools for 
the rest of the paper. In section 3, we present the system that will be studied, and try to 
qualitatively explain the structure of the multiple interface which appears. In section 4, we 
write the equations governing the system and solve them, thus finding the lengthscales and 
concentration profiles which characterize the dynamics of our system (note that, if interested 
only in results, our reader can skip directly to section 4.3). Section 5 closes the article with 
some concluding remarks. 

2 A short overview of two- component interfaces 

We start here by a general overview of the standard theories describing two-species interfaces. 
By no means do we aim here at an exhaustive survey of the theoretical or experimental literature 
on the subject: we simply quickly recall the most basic theoretical elements on which our present 
work builds. 

Let us denote the two polymer species we consider as A and B. We will cover interfaces 
occurring between an A-melt facing a B-melt, as well as between two A-B blends facing each 
other but of differing A/B proportions. We will also mainly put the emphasis on the effect of 
varying the Flory parameter xab between A and B, with the assumption that both A and B 
have the same chain length N. Situations where Aa and N-q differ (e.g. a melt of long chains 
facing a melt of short chains), which are still subject to some debates, will merely be hinted at. 

It is a remarkable point that the kinetics of formation of interfaces between chemically 
different polymers is controlled by thermodynamic as well as kinetic factors. To take into 
account the driving or slowing role played by energetic parameters in diffusion processes, it is 
much convenient to use Flory-Huggins' derivation of the free energy (see eq. 6 for the expression 
of the free energy / in this model, and also note 6): thinking in terms of a lattice model, 
the dimensionless Flory parameter xab characterizes the enthalpy of mixing at the molecular 
level, by comparing interactions between neighboring polymer segments of the same species (i.e., 
association energy eAA for two neighboring A segments, €bb for two B segments) and interactions 
between different species in contact (cab for an A next to a B). If z denotes the coordination 
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number in the lattice, we define 

^(2cab - CAA - £Bb) ,-. s 

XAB = ^ (1) 

In words, xab is the energy change, in units of the thermal energy kT, when a segment of A is 
taken from an environment of pure A and swapped with a segment of B from an environment 
of pure B. The Flory parameter xab can be positive or, much more rarely, negative. If the 
only interactions existing between A and B are van der Waals forces, xab shall be positive. 
Negative values of xab do appear in polymer couples displaying specific chemical interactions, 
such as hydrogen bonds. The three following subsections briefly discuss the different possible 
cases (xab > 0, xab = 0, xab < 0) and their physical consequences. 

2.1 Immiscible components: x AB > 

As stated above, in the absence of specific interactions or structural similarities, the Flory 
parameter is positive and usually ranges from 10~ 3 to 1CT 1 . This is by far the most usual 
situation for polymer pairs. For polymer pairs differing only by isotopic substitution (e.g., when 
deuterium is substituted to hydrogen in the monomer structure), very small, positive values as 
low as xab — 1CT 4 can be found. 7 

In this case of positive XAB) the polymer pairs are immiscible for all but the lowest molecular 
weights: if N = Na = Ab is the length of both polymers A and B, it can be shown that, for a 
Flory parameter xab greater than a critical value xc = 2/ A <C 1, phase separation occurs and 
we end up at equilibrium with macroscopic regions of pure A and pure B, separated by a sharp 
interfacial layer: no macroscopic mixing occurs. 

Let us precise the nature of this interface between the two coexisting phases by a simple 
scaling argument. 8 The characteristic width w of the interface can be estimated by a balance 
between the chain entropy (of order kT per chain), which tends to widen the interface, and 
the unfavorable enthalpy of mixing (of order xab per segment), which tends to narrow the 
interface: thus the typical length Ao of a loop of an A-chain penetrating a B-rich region is given 
by NoXABkT ~ kT, that is, Ao ~ 1/xab- Because the loop has the form of a random statistical 
walk, this contour length corresponds to a (straight) penetration distance w ~ ay/N milx (where 
a is the monomer size). We finally find that the characteristic width of the interface is 

w r (2) 

This simple analysis is in accordance with the results obtained from a more rigorous de- 
scription based on a square-gradient model of the interfacial free energy. 8 Note that the scaling 
analysis leading to this result is valid provided that xab A » 1: equation 2 then predicts the 
mixing region to be much smaller than the polymer coil size, and thus, at the macroscopic 
scale, the interface is to be considered as very sharp one, with a steep profile. This property of 
interfaces between immiscible melts will be used thereafter in our solution (cf. Sec. 3.2). 



4 



2.2 Entropy-driven mixing: Xab = 

We now consider the case of a zero Flory parameter. Rigorously speaking, this situation 
arises when the two polymers put into contact are identical (same chemical structure) . But it is 
also a useful approximation to describe mixing in systems with a very small xab 5 be it positive 
or negative (typically when |xab| *C 1/N). 

Starting from the initial situation where two polymer pieces are put into contact, the final, 
equilibrium state reached by the system is a homogeneous one (the initial interface has totally 
disappeared) with complete mixing of the polymer components. As there is no enthalpy gain 
associated with mixing species (xab = 0), the evolution of the system is driven solely by gains 
in the translational entropy of the polymer chains. 

The dynamics towards this final state can be complex, especially for entangled polymers 
where several regimes of non-Fickean diffusion successively appear. In these entangled polymers, 
the first stages of the interdiffusion process between melts of equal molecular weight have been 
described first within a scaling approach in Ref.10, and then detailed by different authors: 11-14 
the chains are initially segregated on each side of the contact plane, with discontinuous con- 
centration profiles abruptly dropping from unity to zero at the separation. Then, as the chains 
start interpenetrating to form a mixed layer, and because of the specific reptational dynamics of 
the polymer chains, this initial discontinuity survives, with a progressively resorbing amplitude, 
until disappearing when the reptation time T rep is reached. The interface is then said to have 
"healed" completely. This feature stands in obvious contrast with conventional mixing where 
initial discontinuities are immediately smoothed out by the diffusion process. 

For times greater than the reptation time (t > T rep ), the interdiffusion dynamics becomes 
purely Fickean, and is described by the classical diffusion equation 

4> = D s V 2 <j> (3) 

The diffusion coefficient D s , known as the "self-diffusion" coefficient, is independent of concen- 
tration, and for long, entangled chains (A greater than the entanglement threshold A e ), is equal 
to 

As = Ao^fcT (4) 

where Ao is the monomeric mobility. 

The use of eq 3 is valid only for facing melts with chains A and B of equal length (symmetric 
interfaces with Aa = Ae). The case of asymmetric polymer junctions, between chemically 
identical melts of different molecular weights (typically a melt of short chains in contact with 
a melt of much longer chains) proves much more subtle. The evolution of such a system also 
displays several temporal regimes and can essentially be understood as follows: 15 ' 16 the long 
and less mobile chains behave like a gel which is progressively penetrated and swollen by the 
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smaller species; but as the gel of long chains swells, it effectively drags the smaller chains within 
it. The global diffusion profile is thus rather complex and results from the microscopic diffusion 
of the small chains relative to the matrix of longer chains, combined with the global motion of 
this matrix. This global flow of the matrix, known as a "bulk-flow" , 9 is essential for a correct 
description of the system and will reappear in the next section devoted to enthalpy-driven 
species. 

2.3 Enthalpy-driven mixing: Xab < 

Let us now consider the case of mixing between species driven by a gain in enthalpy, 
i.e., when xab < 0. A few dozens of A/B pairs with xab < have been reported, like 
PS/PVME (polystyrene/poly- vinylmethylether) 17 or PVC/PMMA (poly-vinylchloride) /poly- 
methylmethacrylate) . Within the temperature range where xab remains negative, such polymer 
couples are fully miscible for all molecular weights. 

The dynamics of "enthalpy-driven" mixing at interfaces between polymers has been the focus 
of many studies in the past 20 years, and its theoretical description has led to some controversy 
in the case where the system is asymmetric, e.g., when the mobility of the A and B chains are 
different (due to different chain lengths). 

Based upon Onsager's formalism of linear irreversible thermodynamics, 18,19 two main ap- 
proaches were proposed. 

The so-called "slow-mode" theory propounded by Brochard and co-workers 20 ' 21 , assumes 
that the flux of species A is the opposite of the flux of B, i.e., Ja = — Jb- In cases where 
the system is asymmetric, this hypothesis of cancellation of fluxes leads to a mutual diffusion 
coefficient dominated by the mobility of the slowest species — hence the name of the theory. 

On the other hand, Kramer and co-workers 22 explicitly introduce a third species in their 
model, which play the role of vacancies in the system and move about in order to allow the 
motion of A and B chains. Then the fluxes of the A and B components need not be equal 
anymore, and their difference Ja — Jb is compensated by a net flux of vacancies Jy across 
the cell. 23 It is furthermore postulated that the chemical potential of the vacancies is constant 
throughout the system, and this, contrarily to the theory of Brochard et al., leads to a "fast- 
mode" diffusion coefficient, i.e., dominated by the fastest species in asymmetric systems. 

In the recent years, several authors have tried to bridge the gap between these "slow" and 
"fast" theories ; see for instance refs. 24, 25 and references therein. 

Experimental data seem to validate the fast-mode theory. 26-30 Nonetheless, to this day, 
the origin and physical interpretation of the vacancy flux used in the model of Kramer et al. 
or in more recent approaches remains somewhat obscure. Therefore, it is conceptually most 
interesting to note that, in an amended version of her model, 31 Brochard showed that fast-mode 
predictions could be retrieved without resorting to any flux of vacancies. 32 Instead, Brochard's 
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corrected approach shows that, in addition to the individual diffusion fluxes of A and B, there 
exists a global, collective motion involving both A and B; when this is taken into account, 
fast-mode results are found without the need of vacancies. The occurrence of such collective 
"bulk-flows" (as they are called), superimposed to pure diffusion, is indeed known as a common 
feature of asymmetric diffusive systems. 9 In the amended model of Brochard, each chain of 
polymer reptates and diffuses within a matrix composed of all the other chains; physically, a 
bulk-flow occurs because this matrix itself is set into motion, due to the difference of mobility 
between the two diffusing species (in other words, polymer chains reptate in tubes which are 
themselves moving). 

To describe the interdiffusion process of these two different polymers, one usually employs 
the so-called "mutual diffusion" coefficient Dm, which relates the total flux Jj of species i (i =A 
or B) to the gradient of its volume fraction fc: Jj = -DmV^. Let us denote 4> = 4>a = I — 4>b 
the volume fraction of species A. In all the theoretical models mentioned above, the mutual 
diffusion coefficient can then be cast into the generic form: 



Dm = AM-0(1-0)-0 (5) 



where A[(p] is a positive function of (j) (homogeneous to a mobility), which differs from model to 
model, and / is the classical Flory-Huggins free energy per site, 6 given by: 

w = ik log + i^r log (1 " ^ + XAB(p{1 ~ ^ (6) 

In the rest of this article, we will only be concerned with the symmetric case, i.e., with 
polymer components of the same chain length (Aa = Ab = A) and the same monomeric 
mobility Aq. It should be pointed out that, in this case, the conceptual subtleties mentioned 
above do not intervene, and interdiffusion models do agree on the same expression of the diffusion 
coefficient. For long, entangled chains (A > A e , the regime of interest to us), and using the same 
notations as in eq. 4, one has simply A = AoA e /A. From eq 5, the mutual diffusion coefficient 
then reads 20-22 



D M = Ao^fcT -<f>(l 



2XAB 



(7) 



>(1-0)A 

A useful limit in practice is |xab|A » 1: as long as 4> and 1 — 4> are larger than the (small) 
value 4> c = 1/(2|xab|A), Dm can be approximated to 

Dm = D i( j>(l - <f>) with D { = 2|xab|AD s (8) 

(where D s is the entropic, self-diffusion coefficient defined in eq 4). 

Equation 8 shows that the diffusion coefficient in the enthalpy-driven mixing of polymers is 
strongly concentration-dependent. In particular, the vanishing of Dm for <fi = and <f> = 1 leads 
to very unusual diffusion profiles when either of the polymer pieces brought in contact is initially 
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pure in one of the components: When only one side is pure, the mixing width on that side is 
finite, with a sharp edge; if both the initial polymer pieces are pure, the width of the mixing 
region is finite on both sides, and, moreover, the diffusion profile within it is a straight line. 

(Of course, as stated above, the simplified form 8 of Dm does not hold when ^ or 
^ — > 1, and one has in fact to return to the full expression 7 for which Dm does not exactly 
vanish. . . Thus, strictly speaking, the width of the mixing region is never exactly finite, and 
enthalpic diffusion profiles near <p = and <p = 1 cross over to conventional, infinite "diffusion 
tails". But this crossover occurs for </> or 1 — </> of order (f) c = 1/(2|xab|-^V), and in the limit 
where |xab|-W S> 1, 4> c is so small that the diffusion tails contain a completely negligible mass 
of material.) 

With this, we conclude our overview on the different aspects of the theory of two-component 
interdiffusion, and come back to the problem of multi-component interfaces. 

3 Description of the three-component system and qualitative 
approach 

In this section, we present the specific multi-component system that we chose to study, and 
give some qualitative insight into the structure and dynamics that are expected to emerge. The 
next section will deal with these issues in a more rigorous way. 

As explained in the Introduction, the realm of multi-species system and their interfaces is a 
very vast one: our aim here is to explore one "model case" among these possibilities, involving 
three species, and where some interesting features do appear. 

3.1 Description of the system and initial configuration 

We now present the system that we will be concerned with in the rest of this article. The 
initial (t = 0) situation is depicted in Figure 1: the system consists of two pieces of molten 
polymer which have just been put into contact. The surface of contact (the initial interface) 
between the pieces is assumed to be planar, and is located at position x = (the x-axis is drawn 
perpendicular to the surface of contact). It is assumed that the problem is invariant in the two 
remaining spatial directions (parallel to the interface). 

The system contains three polymer components which are initially distributed as follows: on 
the left-hand side, we have a blend of two polymers, A and A*; on the right-hand side, we have 
a pure polymer, i.e., only one species which is called B. The important step (which is really 
defining our system) is how to specify the nature of the pairwise affinities between these three 
components. Our choice is the following: we assume that A and B are attracted to each other 
(through some specific interaction), and are thus miscible in all proportions; on the contrary, 
A* and B strongly repel each other and thus have a tendency to form fully segregated phases; 
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Figure 1: Initial situation of the three-species system. The left-hand side is a blend of two polymers, A 
and A*, with initial volume fractions 4>a — 4>o and (j)A* = 1 — 00- The right-hand side contains only one 
species, called B, with initial volume fraction 4>b = 1. 

and, finally, A and A* are fairly indifferent to each other (i.e., neither attracted or repelled, or 
not much). These features correspond to the following set of Flory parameters: 

Xab <0, with |xAB|A r > 1 (9) 

Xa*b>0, with xa*b^>1 (10) 

XAA*-0, with |xaa*I^«1 (11) 

As we are mainly interested in the effect of the contrast in miscibility between the polymers, 
and not in chain length effects, it has been assumed above that all polymers have the same 
polymerization index: 

N A = N A * = N B = N. (12) 

A three-component system with such features as specified in eqs 9-11 would be quite diffi- 
cult to realize experimentally, as the above requirements seem at first rather contradictory; they 
would certainly be achieved only through a fine tuning of the chemistry of the different poly- 
mers involved. We discuss some practical possibilities in the concluding remarks of the article 
(Section 5). 

Our choice of system is mainly guided here by the fact that it is one case that leads to 
the most interesting interfacial structure and dynamics. Prom the theoretical point of view, 
this interfacial problem is also original in the sense that, as will be seen shortly, it leads to 
simultaneous, coupled, entropy and enthalpy-driven diffusion processes. 

3.2 Qualitative approach to the multiple interface structure 

In this section, we present qualitative arguments explaining why, from the initial situation 
depicted in Figure 1, the interface between the two melts develops a "multiple" structure, i.e. 
displays several layers or regions of different nature. 
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Zone IV: 
A + A* + B 



Zone I: 
A + A* 



Zone II: 
A + B 



Zone III: 
B 



xj(t) 



xil(t) 



X 



Figure 2: Multiple structure of the interface, with several regions, for t > 0. The arrows indicate the 
direction of motion of the regions' boundaries, x = —x\{t) and x — xu{t). The quantities 4>\ and 
indicate the values of the volume fraction in A at the left and right boundaries of region IV. 

Quickly after the onset of contact between the two melts, the interface takes on a config- 
uration as depicted on Figure 2. As can be seen, the interface is divided in four main zones, 
labelled I-IV, and the situation can indeed be intuitively understood as explained now. (We 
do not consider here the transient evolution which takes place in the first instants after contact 
and leads into the configuration of Fig. 2. The interested reader is referred to Appendix A for 
a more detailed discussion of this transient.) 

After the onset of contact, species A and B (which are attracted one to another) start to mix 
together, thus forming a growing mixing layer (region II) around the initial position (x = 0) of 
the surface of contact between the melts. 

On the other hand, A* is repelled by B, and thus, as region II broadens and the A+B 
mixture spreads, A* must recede towards the left. Thus region I, which we define as the domain 
of coexistence of A and A*, is progressively shrinking, with its right boundary moving towards 
the left. 

At the junction between region I (containing A and A*) and region II (containing A and 
B), there must be a region where all three species coexist (region IV). However, we expect this 
region IV to be of very small extension as compared to the others, since we know from section 2.1 
that strongly immiscible polymers put into contact, like here A* and B, tend to form a very 
sharp interface with a steep profile; hence, in the remaining of this article, we will neglect the 
thickness of region IV, and we will not consider any detailed concentration profile within it. It 
must be emphasized however that region IV, albeit tiny, plays an important role in the dynamics 
of the system, since all A chains which want to migrate from region I to II must cross region IV 
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(this will be discussed again in the next section). 

Finally, there is a last region, III, where B is alone: this comes from the fact that, as explained 
in Sec. 2.3, when one of the polymer pieces in contact is initially pure (here, fa = 1 at t = for 
the piece on the right), enthalpy-driven interdiffusion generates profiles with mixing only over a 
finite width; thus, there is on the right side of region II a (time-dependent) location where the 
volume fraction of A reaches zero, and this is what defines the beginning of region III. As the 
mixing region II broadens, the boundary of region III shifts to the right. 

We also introduce in Figure 2 some notations that we will be using intensively in the next 
sections: x = —x\(t) gives the position of region IV, and, in our limit where region IV is taken 
infinitesimally thin, plus or minus exponents (x = — xf) will respectively refer to the right and 
left boundaries of region IV; x = xu{t) gives the position of the boundary between regions II 
and III; and fa(t) and fa(t) respectively denote the volume fractions of A just on the left and 
right of -xi, i.e., fa = 0a| x= _ x - and fa = </>a| x= _ x +- 

After having presented an intuitive description of the interfacial structure and dynamics in 
our system, we will devote the next section to the establishment and solution of the equations 
governing them. 

4 Governing equations and results 

In this section, we first establish the various equations governing the system (Sec. 4.1), and 
then solve them numerically (Sec. 4.2). The physical results of our computations are gathered 
and commented on in Sec. 4.3. 

4.1 Governing equations 

From now on, the volume fraction of A, will be simply denoted by (ft, i.e., we will use 

fax,t) = (f> A (x,i). 

Let us start by enumerating our unknowns: these are the locations of the different regions' 
boundaries (x = —xj and x = xu), the values of the volume fraction in A at x = —xf (fa and 
fa), and the time-dependent profile of fax, t) within regions I and II (by definition, (ft = within 
region III). We will thus need a set of six independent equations to determine the solution of 
the system. 

It should be noted that, from the knowledge of (ft, it is straightforward to retrieve the volume 
fractions of the other components of the system (using the total volume fraction condition 
(ft + 4>a* + fa = 1): in region I, where only A and A* are present, one has <ft A * = 1 — (f> and 
fa = 0; in region 11, fa = I — (ft and (ft A * = 0; finally, in region III, fa = 1 and (ft = 4>A = 0. 

We now establish the six independent equations which govern the dynamics of our system. 
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Diffusion equations 

Our first equations simply describe the diffusive motion of polymer A. We must use two 
different equations, depending on the region considered within the multiple interface. 

In region I, A chains diffuse inside an A+A* mixture, with Xaa* — (eq 11); we are thus in 
a case of entropy-driven diffusion which is accounted for by equation 3, 

<p - D s V 2 <p = (region I) (13) 

and the initial and boundary conditions are as follows: 

4>{t = 0) = <j) ; <j)(x ->■ -oo) = 4> , (f>(x = -Xj ) = 4>i (14) 

In eq 13, the entropic diffusion coefficient D s is a constant parameter (see eq 4). 

In region II, A chains diffuse into an A+B mixture, with xab < and |xab|^ ^ 1 (eq 9); 
we are there in the presence of an enthalpy-driven diffusion, for which the diffusion coefficient 
becomes concentration-dependent (see Sec. 2.3). Recalling eq 8, we have 

D { = 2\ X ab\ND s (Z> f »A,) (15) 
and the diffusion equation in region II can be written as 

<j) - £> f V.(0(l - <A) V 0) = (region II) (16) 
with the initial and boundary conditions 

(j)(t = 0) = 0; 4>(x = -x+) = cf> 2 , 4>{x -► +oo) = (17) 

(see note 33 ). 

Expansion of the A-B mixing region on the right 

Apart from the two diffusion equations 13 and 16, a third equation can be written, which 
governs the expansion rate of the mixing region II (A-B mixture) on the right, that is to say, the 
motion of the boundary x\i(t). By definition, x = x\\{t) denotes, for any time t, the moving point 
where the concentration profile eft reaches zero (and beyond which only polymer B is present). 
As shown in the Appendix B to this article, this very definition implies that the velocity xn of 
this point has a simple relation to the gradient of 4> on its left as follows: 

in = -AV0| - (18) 
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Conservation of the A and B species 

We must naturally also take into account equations that ensure the conservation of the 
different polymer components in the system. 

It is shown in Appendix B that the conservation equation for species A leads to the following 
relation between the velocity x\ of the interface of regions I and II, and the volume fraction <f>\ 
of A and the gradient V(f>\ x= _ x - on the left side of it: 

il = -I3^;- £, - v ^*=-xr ( 19 ) 

The conservation of the B species, on the other hand, gives a relation between the same 
velocity x\, and the volume fraction and gradient on the right side of this interface: 

xi = -<\>i ■ DfVct>\ x= _ x + (20) 

Note that the conservation of the third species, A*, follows automatically from the conser- 
vation of the two others, through the relation between volume fractions (f> + 0a* + <fe = 1, and 
thus does not bring us a new equation. 

The role of region IV 

At this point, we have written down five independent relations; as there are six unknowns, 
there must be a last physical constraint determining the dynamics of the system. We have so far 
totally ignored region IV (as defined on Figure 2), where the three species A, A* and B coexist. 
As stated earlier, this is because this region is expected to be of very limited extension, and, 
as such, has rightfully been discarded, e.g., from the A and B conservation equations discussed 
above. However, despite being negligible in mass, region IV must clearly play an essential role 
regarding the dynamics, since all the A material that diffuses from region I to region II has to 
transit through it. 

Describing precisely the interdependent diffusive process at work within the three-body re- 
gion IV — solving the "inner problem" , in the terms of boundary-layer theory — would undoubt- 
edly represent a complex task and require a study of its own. We rather propose to bypass 
this difficulty, as we in fact only need very limited information for our purpose: we are only 
interested in the effect of region IV on the outer dynamics, not in internal details. 

We make the simple assumption that region IV opposes no particular resistance to the flow 
of polymer A passing through it: region IV is "permeable" to A as other regions in the system 
are, or in other words, the mobility of species A in region IV is not significantly different from 
the mobility in region I or II. This assumption on mobility appears natural if one considers that 
although region IV differs from other regions by significant changes in composition (it is the only 
region where all three species overlap), it does not drastically differ in microscopic structure; 
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hence, the reptational dynamics, which control mobility, should not change significatively within 
region IV as compared to other regions. (The situation would completely be different if region IV 
was, for instance, made of a densely crosslinked network which would hinder the passage of 
polymer chains; one would then rather work under the opposite, "high resistivity" assumption). 

From this permeability assumption, we can deduce a relation on the chemical potential of 
species A on each side of region IV: a chemical potential gradient must be present to drive the 
diffusion of A through region IV, but because this region is both of "low resistivity" and because 
it is very thin (as argued earlier in the text), the overall drop of chemical potential over region IV 
will be minute as compared to the potential drop over other (larger) regions, and can therefore 
be completely neglected. 

Denoting /Ua(— xf ) the chemical potential of species A just on the left of region IV, and 
^a(-x^) the potential on the right, 34 we thus have 



As shown in Appendix B, once the expressions of the chemical potentials are written out, 
this equality between potentials provides our last governing equation in the form of a relation 
between the volume fractions 4>\ and 4>2 at the borders of region IV: 



We are now in possession of six independent equations which govern the dynamics in our 
system, and can proceed to solve them. 

4.2 Technical solution 

We now present the procedure followed to solve the set of governing equations derived in 
the previous section (eqs 13, 16, 18, 20, 19, 22). We will here focus on the technical aspects of 
the solution, while, for the sake of readability, we have grouped together the presentation of the 
results and their physical content in the next section (Sec. 4. 3). 

The major difficulty attached to the interdiffusion problem as defined by this set of equations 
is that it involves two moving boundaries, x = —x\(t) and x = x\i{t): the solution of the diffusion 
equations in the system requires the application of boundary conditions defined at the moving 
boundaries, but the motion of the boundaries themselves are directly determined by the diffusive 
solutions. This self-consistent nature of moving boundary problems make them generally much 
less straightforward to solve than a conventional diffusion problem; we will here follow the 
solution scheme described in Crank. 9 

Our problem falls into the so-called "Class A" category, 9 characterized by the fact that 
the motion of the moving boundaries is uniquely due to the transfer of diffusing substances 
across them. (This is indeed the meaning of the conservation equations 38 and 40 discussed in 



Ma(-^i ) = VA(-xf) 



(21) 



01 = <t> 2 exp[-|xA B |A r (l-0 2 ) 2 ] 



(22) 
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Appendix B). It is then possible to show that, in a generic fashion, the solution will have the 
following properties: 9 

(i) The concentration of the diffusing species at the moving boundaries is constant, i.e., 4>\ 
and 4>2 are independent of time. 

(ii) The motion of the moving boundaries is diffusive, i.e. x\(t) ~ {D s t) l l 2 and x\i(t) ~ 
(At) 1/2 - 

We build the solution to our problem with the help of these useful results. In accordance 
with property (ii) above, we introduce the new quantities a and b defined by: 

xi(t) = 2a{D s t) l/2 (23) 
x u (t) = 2b{D i t) 1 ' 2 (24) 

The numerical value of a and b will have to be determined as a result of our solution. We also 
define 

e =S=ra «> < 25 > 

which is a small parameter (see eqs 8 and 9). 

In the limit where e <C 1, equation 22 simplifies considerably: 

* = *^(-^r £ )^° < 26 > 

Thus, as long as (1 — 02) S> e 1 / 2 (i.e., as long as <p2 is not too close to unity), we can as well 
consider that (p\ = for the rest of our solution to a very good approximation. 

We can then immediately solve the diffusion equation 13, valid in region I, along with the 
boundary conditions 14. The solution can as usual be looked for as a function (f>(x,t) = f(u) of 
the reduced variable u = x/2\/D s t, and from eqs 23 and 26, we see that the boundary condition 
on the moving boundary simply becomes <p{— xj~,i) = /(— a) = 0. Solving the equation on /, 
we finally find the following concentration profile in region I, which involves the error function 
erf u = (2/- v /7r) J q m e~" 2 du, classical for semi- infinite media: 

x 



4>{x, t) = 4> ' 



1 — erf a 



1 + erf 



2JlXi 



(region I) (27) 



We next turn to the determination of a, and thus to the determination of the position of the 
moving boundary at x = — x\. Equation 19 will serve this purpose: with <p\ = 0, it simplifies to 
x\ = — D s X7(f)\ x= _ x - . Then, with the help of eq 23, which yields x\ = a(D s /t) 1/ ' 2 , and of eq 27 
which yields an expression of V<p\ x= __ x - , eq 19 finally brings the relation 

v / ^ae a2 (l-erfa)=0 o (28) 

This equation on a can be solved numerically to give the value of a as a function of the fundamen- 
tal parameter 4>q (which describes the initial state of the system) . The result of the numerical 
solution is shown in Figure 7. 
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The next step is to solve the diffusion equation 16 under the conditions 17, in order to find 
the concentration profile in region II. As equation 16 is a nonlinear equation without analytical 
solution (to the best of our knowledge), we will resort to numerical integration. We here again 
look for a solution 4>(x,t) = g(v) of the (new) reduced variable v = x/2^D{t. The resulting 
differential equation on g turns out as: 

and the boundary conditions of eq 17 translate to 

(i) (f>\ x= - x + = g\v=~a^i = <h, (h) 0U-+OO = Slv-H-oo = ( 30 ) 

It has to be noted that these boundary conditions, as such, are not sufficient to directly proceed 
to numerical integration, because the quantity <p2 appearing there is also an unknown. We thus 
need a supplementary condition, which will be provided by the conservation equation of the B 
species (eq 20): using the definition of a (eq 23) and the definition of the reduced variable v, 
eq 20 can be rewritten as 



dg 

dv 



2a ^ € (31) 



v=-a^fe (p2 

We have now enough equations to compute g(v) and the associated concentration profile 
in region II. We solve the set of equations 29-31 by a trial-and-error procedure: we choose 
an arbitrary value for the unknown 4>2 ; based on the corresponding initial conditions provided 
by eqs 30- (i) and 31, a numerical integration of eq 29 is performed; the obtained solution g is 
then evaluated at infinity, and its value compared to zero (eq 30- (ii)); if different, the whole 
numerical procedure is resumed with a new chosen value for 02, until the correct value of 4>2 , 
satisfying eq 30-(ii), is found. At the end of this process, we eventually obtain the actual profile 
of concentration cj)(x, t) = g(v) in region II, along with the value of <p2 in our system. (See 
further down for results.) 

We are at this point in possession of the profiles (p(x, t) in regions I, II and III, as well as 
of the numerical values of a, <p\, and <f>2- To complete our solution, we still need to find the 
value of b, which is related to the position of the border between region II and III. This is easily 
performed numerically: by transposing the equation of motion for xn established earlier (eq 18) 
in terms of the self-similar function g(v), and applying the definition of b (eq 24), one obtains 
the relation 

dg 



dv 



(32) 

v=o 



i.e., b is a quantity such that the slope of g at v = b is equal to twice the value of b (with 
changed sign). The search for such a b can be carried on numerically with the function g 
calculated previously: the graph of b as a function of 4>o can be found on Figure 7. 
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Figure 3: Sketch of the generic self-similar diffusion profile for species A: the volume fraction <\> has been 
drawn versus the reduced variable u = x/2y/D s t. The volume fractions used in the text, (f>o and (fa, 
are also represented. In domain I (i.e. for u < —a), the concentration profile is an error function; in 
domain II (i.e. for —a<u< 6/e 1 / 2 ), the profile is quasi-linear. In region III, A is absent. 

4.3 Presentation of the results 

We now present and comment on the results of the model, as obtained through the procedure 
described in the previous section. 

As seen in eqs 23 and 24, the different "compartments" or layers composing the multiple in- 
terface are separated by diffusive fronts: x = —x\ (t) = -2a(D s t) l l 2 and x = x\i(t) = 26(Dfi) 1 / 2 , 
with a and b are numerical parameters (see below for their values). The evolution in time of the 
concentration profiles of the different components is thus self-similar, and for this reason, it is 
useful to present them here in terms of the reduced variable u = x/2(D s t) 1 / 2 . 

We remind our reader that there are two main physical parameters to our solution: <fo, 
which represents the content in A of the initial A+A* blend, and the small parameter e = 
1/(2|xab|A0 *C 1, which is related to the strength of the attraction between A and B chains 
and the polymer chain length. 

As these parameters were varied over a wide range of values (4>q between 0.1 and 0.9, e 
between 10 -5 and 0.1), it was observed that, remarkably enough, the (self-similar) diffusion 
profiles obtained retain the same aspect throughout; the main effect of varying the parameters 
is simply to rescale one or the other of the characteristic dimensions of the self-similar solution 
(i.e., stretch or shrink it in some region). 

This generic self-similar diffusion profile for the A species is sketched in Figure 3 in terms 
of the reduced variable u. The fronts x = —xj(t) and x = x\\(t) are respectively located at 
u = —a and u = b/y/e. Starting from the left in the Figure, the profile in region I is given by an 
"error function", typical of diffusion problems, and drops to zero at the approach of region II 
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Figure 4: Numerical profiles for the volume fraction (f>(u) of the A species (with u = x/2\JD s t). The 
chain length and Flory parameter are N = 10 3 and xab = 0.1, yielding e = 5 • 10~ 3 . The solution (f>(u) 
has been represented for two values of the initial volume fraction of the A species in the A-A* blend: 
4> a = 0.3 and O = 0.7. 

(since we proved that 4>\ ~ in most practical cases) . At the boundary between regions I and II 
(we neglect the thickness of the intermediate region IV), there is a discontinuous concentration 
jump from (f> = to <j> = fo, and then the diffusion profile decreases smoothly until reaching 
zero again at the border with region III. It should be emphasized as a very interesting feature 
of the solution that the amplitude 4>2 of this concentration jump remains constant in time and 
does not die away during the diffusion process. 

Another interesting feature is the depletion near u = — a of species A (0 — > 0). The depletion 
is a direct reflection of the fact that it is energetically much more favorable for A to be in the 
A+B blend of region II (xab < 0) than in the A+A* blend of region I (xaa* — 0); therefore, 
as time passes, region II continuously sucks A-chains out of region I, creating in the latter a 
depleted zone which grows with time (as \JD s t). 

We also note that, from an empirical point of view, the profile in region II does not have a 
very significant curvature; if needed, we may thus approximately regard it as a straight line. 

As an illustration of the generic sketch of the diffusion profile described above, we have 
plotted in Figure 4 two actual profiles, computed numerically for different values of </>o (and 
same e), which do show that, within some rescaling, the generic shape of Figure 3 is retained in 
both cases. 

For completeness, we also give the corresponding generic self-similar profiles for the other 
mixture components in Figures 5 and 6. As reminded in Figure 3, the diffusion profile involves 
several characteristic quantities a, b, and 02, whose values we now consider in detail. 
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Figure 5: Sketch of the generic self-similar diffusion profile for the A* species: the volume fraction (f)p* 
has been drawn versus the reduced variable u = x/2^/D s t. In region I, (j) A * = 1 — <fi and is thus given by 
a complementary error function. In the other regions, A* is absent. 




Figure 6: Sketch of the generic self-similar diffusion profile for the B species: the volume fraction 0b has 
been drawn versus the reduced variable u = x/2y/D s t. In region I, B is absent; in region II, the profile 
of B is the complement of A (0 B +0 = 1) an d is quasi-linear; in region III, B is alone. 
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Figure 7: Values of the numerical factors a and b, which characterize the two moving boundaries x\{t) = 
—2a(D s t) 1 / 2 and x\\{t) = 2b{Dft) 1 / 2 , plotted versus the initial volume fraction (f> a . (cf> a is ranging from 
0.1 to 0.9 and e is fixed at 5 • 10~ 3 .) 

In Figures 7 and 8, we plot the variations of these quantities as functions of our parameter 
4>o, i.e., a((po), b((fro), and 02(0o)> at fixed e (i.e., we work at given system chemistry, and change 
the initial composition of the system). It is observed that these quantities show a steady increase 
with 4>o. Results are presented for the range 0.1 < cfro < 0.9, which should cover most practical 
cases; see note 35 about situations where cfro — > or <fro — > 1. 

It is also possible to give rough analytical estimates of the dependence of a, b, and 4>2 with 
respect to our other main parameter, e. The quantity a is given by the solution of eq. 28, which 
involves only (fro; thus we see that a has no functional dependence on e (at least in the regime 
e <C 1 considered in this article). Let us now estimate the dependencies of <fr>2 and b on e: for 
that purpose, we will consider, in first approximation, that the concentration profile is linear in 
region II (see above). Writing as previously the solution (fr(x,t) = g(v) in terms of the reduced 
variable v = x/2(Dft) 1 / 2 , the ("constant") value of the slope dg/dv is easily estimated: since 
g(v) goes from 4>2 at the left boundary of region II (v = —ay/e) to zero at the right boundary 
(v = b), we have dg/dv ~ — fo/ib + ay/e) ~ — fo/b. We then use this expression of the slope into 
eq 30- (i) and eq 31, and solve this set of two equations for b and (fri- This brings the following 
estimates: 

/a 2 \ 1/6 

b*(je\ , 2 ^(2a 2 e) 1/3 (33) 
In terms of scaling laws with respect to e, we thus have 

b ~ e 1 / 6 , 02 ~ e 1/3 , a ~ e° (34) 

(where the last equation is meant to recall that a has no dependence on e). 

Finally, we conclude the presentation of our results, by emphasizing one unusual and notable 
feature of our solution: as time goes on, region II, where A and B mix, grows in a very asymmetric 
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Figure 8: Plot of the quantity 02, which corresponds to the amplitude of the concentration jump at the 
border between region I and II, versus the parameter <j)Q. (<fio is ranging from 0.1 to 0.9 and e is fixed at 



fashion, relatively to the position of the initial interface (x = at t = 0); this is due to the fact 
that its two boundaries are of very different nature, with the right boundary x = x\\ having a 
much faster diffusive motion than the left boundary x = —x\ (because D$ » D s ). The ratio of 
these two boundaries' position is in fact independent of time and can be estimated as 



The sluggishness of the motion of x\ relatively to x\\ thus becomes more pronounced as e becomes 
smaller, i.e., if the attraction between A and B is made stronger or if the polymer chains length 
is increased (see eq 25). 

We have now completed the solution of the dynamics of our three-species system, and have 
presented the results that were obtained. 

5 Concluding remarks 

In this article, we studied the evolution of the interface between two polymer melts, in a 
specific case where three species of strongly contrasting chemical affinities were involved. It was 
found that, due to the simultaneous presence of entropy and enthalpy-driven diffusion processes, 
the dynamics at the interface is unusual (with a very asymmetric growth of the mixing layer), 
and that moreover, the interface shows a peculiar spatial structuration in three distinct layers 
of different chemical compositions. 

Let us now close with some remarks. Obtaining a triad of polymers with chemical properties 
such as specified in eqs 9-11 is certainly a difficult task. One potential way could be the following: 



5- 10- 3 .) 




(35) 
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(i) having A and B strongly attracted to each other through the existence of hydrogen bonds, 

(ii) having A* a very slightly modified variant of A, where hydrogen bonding to B becomes 
unavailable; then between A and B, van der Waals attractions will dominate, leading to repulsion, 
while hopefully A and A* will remain similar enough and roughly indifferent to each other. 
Another, powerful way of achieving such a system may also be through the use of copolymers. 

The specific case which has been described and studied in this article is, as has been argued 
in the Introduction, only one limiting situation in a vast range of multiple species interfaces 
in polymer systems, and has been selected for the strong opposing tendencies in the chemical 
affinities of the involved polymers. Our hope is however to demonstrate on this example the 
richness of multiple-species interfaces, where many other such "selected cases" would be worth 
investigating and studying for their peculiar structure and dynamics. 

A related, though different, problem to that of two multi-component melts facing each other, 
is the situation of a melt facing a loosely crosslinked gel (this perspective was suggested to us 
by one of the referees). In the example studied in this article, an A+A* blend would be brought 
into contact with a gel of B. The subsequent evolution might develop interesting features: the 
motion of the x\\(t) front (as defined on Fig. 2) would describe the imbibition of the B gel by 
the A species, while the displacement of the xj(t) border would be dictated by the swelling 
of the gel — this would imply the inclusion of new elastic terms in the model's equations. We 
note however that, in practice, the dynamics of polymer diffusion into a network can prove very 
subtle, with a predominant role of preexisting heterogeneities (e.g., in crosslink density) within 
the network. 36 

From the application point of view, the ability to spontaneously form (at least in some cases) a 
multi-layered interface between two polymer pieces, might lead to interesting perspectives in the 
design of multilayer polymer films 37 — which are widely used for product packaging, for instance. 
Whereas conventional technology includes processes like coating, lamination, or coextrusion to 
form multilayered films, here the layering appears by means of a purely diffusional process. The 
chemical composition and the spatial organization of the layering could be controlled through 
an appropriate choice of chemical affinities between the different species in interaction, and once 
the desired result is attained (e.g., prescribed layer thicknesses), the system could be quenched 
to freeze the diffusion process. 
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Appendix A. Transient behavior at the onset of contact 



In this Appendix, we discuss the transient behavior taking place in the first instants after 
the contact between the initial melts is established, and how it leads to the layered configuration 
depicted on Figure 2 of the main text. 

Let us consider the situation immediately after the onset of contact (see Fig. 1), from the 
point of view of species B: due to thermal agitation, a few B molecules start to extend chain 
portions into the A+A* blend, of respective volume fractions 4>a and 0a*. These first B chain 
portions will then experience an A+A* environment with an average, "effective" Flory parameter 

Xb - 0AXAB + 4>a*Xa*b (36) 

(in a mean-field picture). Given the opposite signs of xab and Xa*B) this effective Xb m ay be 
positive or negative depending on the composition of the A+A* blend. Defining 

4>c = | Xa*b I / (xab + | Xa*b I ) (37) 

we have two cases depending on the initial (t = 0) content <px = <fo in A of the A+A* blend: 
either 4>o < 4> c or 0o > 4> c . We discuss below the two scenarios that ensue for the transient 
evolution of the system at short times. 

If the initial situation is such that </)q < <f> c , we see from eq 36 that x B > 0, which means 
that the environment seen by the B chains exploring the initial A+A* blend is enthalpically 
unfavorable. Therefore the triple region A+A*+B created by this diffusion remains of very 
small extension and corresponds to our region IV in Fig. 2. On the other hand, the diffusion of 
species A into the B melt is very favorable (xab < 0), and thus A molecules do at the same time 
migrate from the initial A+A* blend through the triple region IV, and create a mixing layer 
A+B as seen in Fig. 2. At that point in time, the layered structure of the interface is established 
and the subsequent dynamics of the system crosses over to that described in the main text (see 
Sec. 3.2). 

Alternatively, if the initial situation of the system is such that 4>o > 4> c , eq 36 tells us that 
Xb < and in that case, the environment seen by the B molecules initially in contact with the 
A+A* blend is favorable. Therefore, it is expected that after the onset of contact, B molecules 
will start diffusing into the A+A* blend and thus a growing region where the three species A, 
A* and B overlap will form. This may at first seem at odds with the layered configuration 
of Fig. 2 and especially with the assumption of a thin triple region invoked in the main text. 
However, if one considers the point of view of species A and A*, simple energetic arguments 
make it clear that the initial growth of the triple region will be short-lived: firstly, the triple 
region makes an unfavorable environment for A* (the local effective Flory parameter for A*, as 
an average of Xaa* — and Xa*b > 0> is always positive) and thus A* will demix and migrate 
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back to the A+A* melt on the left; secondly, rather than staying within the triple region where 
part of their contacts are made with A* molecules (with no enthalpy gain), A chains will seek 
to augment favorable xab < contacts and will thus preferably migrate into the pure B region 
on the right of the triple region — thereby creating the A+B region described in the main text as 
region II. As both these processes occur at the expense of the triple region, the size of the latter 
can only remain modest. Qualitatively at least, it thus becomes apparent, that here again, the 
system rapidly evolves to the configuration of Fig. 2. A more precise description of the transient 
dynamics just sketched would nevertheless be desirable and remains ahead of us. 

(Note: in the limit where the initial A+A* blend is in fact almost purely composed of A, 
i.e., 4>o very close to unity, we expect some of the above arguments to break down; the system 
will not form layers at the interface, but should simply evolve as in the known situation of a 
pure melt of A in contact with a pure melt of B. 20 ) 



Appendix B. Derivation of governing equations 

In this Appendix, we give a detailed derivation of some of the governing equations presented 
in section 4.1. 

Equation of motion for Xu(t) 

Equation 18 is obtained as follows. By definition, xu(t) is the point such that, at any time 
t, (f)(xu(t),t) = 0. Since is a constant over time at that point, we also have that the total time 
derivative at that (moving) point, d(p/dt\ x=Xu , must be zero: 



~dt 



= <j>\x=x u + X U .V(j)\ x=Xu = 

x=xu 



Furthermore, using the diffusion equation 16 and <p\ x=Xu = 0, another expression for (ft is easily 
obtained: 

<j>\x=x u = DfC^(p\x=x u ) 2 
Substitution of the latter equation into the former yields equation 18. 

Conservation of the B species 

We now show how relation 20 results from the conservation of the B species. 
As a little reflection from the consideration of Figures 1 and 2 can make clear, the conserva- 
tion of the B species reads as follows: 

/XII (t) 
<j) B dx = x ll (t)-l (38) 
-xi (t) 



24 



This expression simply states that, when considering the situation at a certain time t > 0, all 
the initial B material which was comprised between x = and the present position x = x\\(t) (a 
quantity equal to x\\ ■ 1 since the initial volume fraction of B was one) has been redistributed by 
virtue of the interdiffusion process all over region II (hence the term 4>b dx on the left-hand 
side). Using <\>b = 1 — <f>, one then obtains 

rxn(t) 

>dx = xi(t) (39) 



-XI (t) 

Differentiating eq 39 once with respect to time, using eq 16 and <p(xu) = 0, one obtains a local 
version of this integral equation: x\ = — </>2 • DfV(f>\ _ +, which is exactly equation 20. 



Conservation of the A species 



We now derive equation 19 from the conservation of species A. The conservation equation 
for A has the form 

/xn(t) r-xn(t) 
<j)&x = (f H) x l {t)+ / (fo-^dx (40) 

-Xl(t) J— CO 

This expression has the following meaning: at time t, the A material enclosed in region II 
(Sx=-xi & comes from the material which was initially present between x = —x\ and x = 
at volume fraction <j>o, plus the amount that has diffused from region I (/^"(^o — 4 1 ) dx). Using 
eq 39 into eq 40, then differentiating with time and substituting with the diffusion eq 13, one 
finally obtains the following local version for the conservation of the A species: x\ = — (1 — 
(pi)" 1 ■ D S V(/>\ X= _ X -, which is precisely equation 19. 

Relation between <pi and <p2 

We now explain how one obtains equation 22 from 21. 

Using classical formulae 1 ' 5 ' 38 to compute the chemical potential of A, one obtains the ex- 
pression 

^ A = fcT[log0 + X AT(l-0) 2 ] (41) 

This expression for holds for blends made of two polymer species, characterized by a Flory 
parameter x (and under the assumption that these species have the same chain length N). Thus, 
eq 41 should not be used to compute within region IV, where three component coexist; but 
it can indeed be used at the border with region I on the left, and at the border with region II 
on the right — there one of the components drops to zero concentration, thereby leaving us with 
two species only (see also note 34 ). 

On the left border of region IV (x = — xj~), we only have the two species A and A* (since 
— ^ 0), and accordingly, we can apply formula 41 with \ = Xaa* — 0; which yields 

A*A(-zf) 



kT 

25 



log 01 



Similarly, on the right border (x = —xf), only A and B coexist, and we find 

Writing the equality between the above two expressions of the chemical potential (as required 
by eq 21) yields the governing equation 22. 
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